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Abstract: We revisit the CDF data on di-muon production to impose constraints on 
a large class of Z' bosons occurring in a variety of Eq GUT based models. We analyze 
the dependence of these limits on various factors contributing to the production cross- 
section, showing that currently systematic and theoretical uncertainties play a relatively 
minor role. Driven by this observation, we emphasize the use of the Bayesian statistical 
method, which allows us to straightforwardly (i) vary the gauge coupling strength, g', of 
the underlying (ii) include interference effects with the Z' amplitude (which are 

especially important for large g'); (iii) smoothly vary the charges; (iv) combine 

these data with the electroweak precision constraints as well as with other observables 
obtained from colliders such as LEP 2 and the LHC; and (v) find preferred regions in 
parameter space once an excess is seen. We adopt this method as a complementary 
approach for a couple of sample models and find limits on the Z' mass, generally 
differing by only a few percent from the corresponding CDF ones when we follow their 
approach. Another general result is that the interference effects are quite relevant if 
one aims at discriminating between models. Finally, the Bayesian approach frees us of 
any ad hoc assumptions about the number of events needed to constitute a signal or 
exclusion limit for various actual and hypothetical reference energies and luminosities 
at the Tevatron and the LHC. 
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1. Introduction 

The search for new physics beyond the Standard Model (SM) is one of the main ob- 
jectives of the current and future collider experiments. A promising signature of such 
physics are neutral Z' bosons which appear in numerous models extending the SM 
gauge symmetry group by an additional 11(1)' factor (for reviews, see [1, 2, 3, 4]). This 
is not only interesting in its own right. In fact, many (or even most) theories, scenarios, 
and models beyond the SM have been supplemented by extra 11(1)' symmetries, in or- 
der to cure or ease specific problems that have arisen there. Thus, by finding a Z' and 
reconstructing the underlying t/(l)' charges, one may obtain clues regarding the un- 
derlying physics and principles, like supersymmetry, large or warped extra dimensions, 
other strong dynamics, etc. 

There is currently no experimental evidence for a Z', not in the electroweak preci- 
sion data (EWPD) [5, 6], nor in the study of interference effects at LEP 2 [7], nor in 
searches for resonance production at the Tevatron [8] or the LHC [9, 10]. It is therefore 
customary to set lower limits on the Z' boson masses, M^/, for a number of models 
and relative to a fixed value of the U{iy gauge coupling, g'. 

However, one can extract more information from all the available experimental re- 
sults than is reflected in a collection of mass limits. Indeed, the information is based on 
a variety of different channels and observables [e.g., cross-sections and asymmetries) 
which can be used to disentangle the underlying model parameters and to diagnose 
the Z'. Moreover, precision constraints come from the Z-pole and related observables, 
from low energy measurements, and from the flavor sector. Thus, one should find a 
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framework allowing to discuss these very different sources simultaneously and trans- 
parently. Bayesian data analysis, being particularly suited for parameter estimation 
(as opposed to hypothesis and model testing), proves very convenient to achieve this 
goal, and in this paper we take a first step in this direction. Specifically, we consider 
as an example the recent di-muon results by the CDF Collaboration [f 1]. 

Working towards the above mentioned goal, we exploit here the implicit features of 
the Bayesian approach i) to study the effects of interference of the Z' boson with the 7 
and Z bosons; and ii) to project exclusion limits for current, future, and hypothetical 
colliders and various luminosity reaches. For the latter, this approach avoids any ad 
hoc assumptions about how many observed or expected events would constitute an 
exclusion or a discovery. 

In our previous study [5] of Z' bosons we analyzed the most recent EWPD, which 
— as usual — was based on least-x^ fits, i.e., the likelihood is given by £ewpd = 
exp(— Further factors Li entering the total (posterior) probability density need 
to be constructed for each data set, 

Pposterior(M^', 5'', ^ZZ', •••) = Ppr\or{M z' , q' , QzZ', •••) X £eWPD X J]^ £i, (1.1) 

i 

where Qzz' is the mixing angle between the Z' and the ordinary Z boson and Pprior is a 
non-informative prior density. We usually take it to be flat for variables defined over 7?." 
having in mind an infinitely wide multivariant Gaussian distribution. In other 
parameter transformation may be in order. The dots refer to further Z' parameters such 
as those characterizing the U (1)' charges. We believe that computing the full ^posterior is 
a (long-term) task worth its effort to make full use of all experimental (and theoretical) 
information. This is particularly obvious when a signal is seen in one or several places, 
and one wants to narrow down the space of possible underlying t/(l)' symmetries, as 
mentioned above. Notice that ^posterior can easily be updated by regarding it as a (now 
informative) prior density and multiplying it with a new factor of Li (we are assuming 
here that there are no experimental or theoretical correlations between various factors 
Li so that the factorization property holds). 

Since among the final goals of our approach is the ability to discriminate among 
models (charges) we review in Section 2 a large and popular class of Z' models based on 
the gauge group and try to put these on an equal footing. Z' models with the same 
charges (at least as far as the SM fermions are concerned) also arise from a bottom-up 
approach [12] when demanding the cancellation of gauge and mixed gauge-gravitational 
anomalies in supersymmetric extensions of the SM together with a set of fairly general 
requirements such as allowing the SM Yukawa couplings, gauge coupling unification, 
a solution [13, 14] to the //-problem [15], the absence of dimension 4 proton decay as 
well as fractional electric charges, and chirality (to protect all fields from acquiring very 
large masses). 
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In Section 3 we lay out the theoretical framework for Z' hadro-production, where 
some technical details relevant to the calculation of cross-sections are given for com- 
pleteness in Appendix A. As a check, we reproduce the CDF limits following their 
approach and extend the limits to models not considered in their original analysis. 
Moreover, we project Mz' limits for various integrated luminosities and center-of-mass 
(CM) energies for the LHC. 

Finally, in Section 4 we discuss the proposed Bayesian statistical method. We derive 
mass limits and compare them with those in Section 3. We then compute exclusion 
contours for some illustrative models, emphasizing the role of interference between the 
Z' and SM amplitudes, especially for large g' . We conclude in Section 5. 



2. The model class 

As mentioned in the introduction, a very large class of f/(l)' symmetries underlying 
the Z' bosons are subgroups of -Ee, and can be written in the form, 

ry, V I • nry , . nry Cl Z R + y/3 Z + C3 Z L^) 

Z = cosa cosp + sm a COS p Zy + smp = ^ . (2.1) 

Here — 7r/2 < /3 < 7t/2 is the mixing angle between the U{1)^ and U{1)^ maximal 
subgroups defined by [16] 50(10) ^ SU{5) x f/(l)^ and Eq -> 50(10) x f/(l)^, 
respectively, and — 7r/2 < a < ir /2 is non-vanishing when there is a mixing term [17] 
between the hypercharge, Y, and U{1)' field strengths oc F^'^F'^^, and this kinetic 
mixing term has been undone by field redefinitions. The f/(l)y, f/(l)^, and U{1)^ 
groups are mutually orthogonal in the sense that Tr(QiQj) = when the trace is 
taken over a complete representation of Eq. The second form appearing in Eq. (2.1) 
uses a different orthogonal basis, U{1)r, U{1)r-^, and U{1)li, which are the maximal 
subgroups [16] defined by SU{3)l,r ^ SU{2)l,r x U{1)l,,r, and SU{2)r U{1)r, 
referring here to the trinification subgroup [18] of Eq — )■ SU{3)c x SU{3)l x SU{3)r. 

The U{iy charges of the particles appearing in the fundamental representation of 
Eq are shown in Table 1 in terms of the parameters ci, C2 and C3, satisfying, 

, C1 + C2 + C3 sgn[| ci - (C2 + C3)] 

tana = — ^= ^= , tanp = — (C3 — 02). (2.2) 

Vi^l- Vi(c2 + C3) yic? + (C2 + C3)2 

The values of a, (3, and the q for some specific models are given in Table 2. We also 
display the charges for the models listed in Table 2 more explicitly in Table 3. The 
general classification of all models with integer charges, as well as models arising from 
breaking chains involving maximal subgroups is the subject of Ref. [19]. 

There are also classes of models described by one continuous parameter. For exam- 
ple, one can restrict oneself to f/(l)' subgroups of 50(10), i.e., those perpendicular to 
the U{1)^ and therefore with C2 = C3 <^=^ (3 = 0. These models are equivalent (up to 
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Table 1: Charge assignment for the left-handed multiplets contained in a 27 dimensional 
representation of Eq. The upper part of the table corresponds to the 16 dimensional rep- 
resentation of SO{10), while the lower part shows the 10 (with an extra anti-quark weak 
singlet, D, of electric charge —1/3 and an additional weak doublet, L, as well as their SM- 
mirror partners) and the 1 (a SM singlet, S). This represents one fermion generation, and 
we assume family universality throughout. The correct normalization (i.e., the one which is 
directly comparable to the usual normalization of the gauge couplings of SU{3)c and SU(2)l 
of the SM) of these charges is obtained upon division by + 3 (c^ + C3). 



non-chiral sets of fermions) to the ones described by the real parameter x and denoted 
by U{l)q+xu in Ref. [30], i.e., with charges defined by Q'^ = —xQ'g = xQ\/3, when one 
identifies, 

[3 x + 1 , , 

tana = W -. (2.3) 

\ 2 x-A ^ ^ 

This class [31] contains the Z' models based on left- right symmetry, Z^ji, which can be 
seen from the breaking, 50(10) ^ SU{3)c x SU{2)l x SU{2)n x U{1)b-l- Inciden- 
tally, Table 2 shows that U{1)b-l is the diagonal combination of U{1)li and U{1)r-^ 
(and indeed was initially dubbed U{1)l+r [20]) and has manifestly left-right symmetric 
charges. However, left-right symmetry is broken in the SM. In the fully 5*0(10) sym- 
metric case, i.e., when all gauge couplings are equal, g = Ql = 9r = 9b-l = Qy (the 
tilde denotes 50(10) normalization which we use here to simplify the discussion), the 

Zlr = cos 9lr {-Zb^l) + ^inOiRZn (2.4) 

must be orthogonal to the Zy since these are obtained by an 50(2) rotation of Zb-l -L 
Zr. This yields the Z^ and tan^^ij = ^2/3 ~ 39°. At lower energies, renormalization 
group (RG) effects will generally split the gauge couplings. One then has [3, 32] 

tan^,« = \/|^ = \l\^4-^ = M^otHw-l. (2.5) 
V 3 qb-l 9y \ 9l 

where the second step uses the relation [31] 5^y^ = 2^^^^ + 3g~^, and where the 
weak mixing angle 9y/ = arctan((7y/(7/,) appears in the last. Thus, assuming manifest 
left-right symmetry, gi = Qr, one finds 6lr ~ 57°. Formally one has the entire range. 
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Table 2: The values of the Cj and (a, /3) parameters for various motivated Z' bosons, 
most of them appearing in the literature (if referenced) . The and the Z^ are bosons which 
do not couple — at vanishing momentum transfer and at the tree level — to protons and 
neutrons, respectively. Similarly, the Z^, Z/, and Z^ bosons are blind, respectively, to SM 
leptons, up-type quarks, and down-type quarks. The Zb-l couples purely vector-like while 
the Z^ has only axial-vector couplings to the ordinary fermions. The overall sign of the Cj 
for each model is not physical, and can be absorbed into the definition of the U{iy gauge 
field when allowing both signs for the mixing angle, Ozz'- However, the sign convention of 
the charges becomes significant once a sign convention for 9zz' has been adopted. Where 
applicable we follow the sign conventions of Refs. [3, 5] and otherwise the more systematic 
sign convention detailed in Ref. [19]. 

^ ^ Glr < 90°, but reahstic breaking patterns [32] suggest 35° ^ Olr 42°. Finally, 
this range can be extended to include all 5'O(10) models by identifying 9lr = a + 
arctan a/2/3. 

Similarly, there is a class of models perpendicular to the U{1)y and therefore with 
Cl + C2 + C3 = <^==^ a = 0. Under the identification, 

tan^^/l^, (2.6) 

these models are equivalent to those denoted by U (l)io+x5 [30], i.e., with charges related 
by <5j = xQ'j^ = xQ'g. Finally, Ref. [30] discussed another one-parameter subset of 
models, U{l)ci-xu, which can be obtained by demanding Q'^ = <^=^ C3 = and 
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Table 3: Explicit U (1)' charges for the models defined in Table 2. The normalization imposed 
by E% symmetry is also shown, and we will use this normalization throughout in order to 
avoid spurious factors in the discussion in the text. Columns within single lines fill out SU (5) 
multiplets, while double lines enclose full SO(ASS) multiplets. 



Q'^ = —xQ'g, and by identifying, 

2^6 X \/3 (x - 1) sgn(x - 5) 

tana = , tanp = . 2.7 

x-5' V5x2 - 2x + 5 

Of course, any other one-parameter subset of models may be considered. All models 
are guaranteed to be free of anomalies due to the absence of an independent cubic 
Casimir invariant from Eq. On the other hand, the U{1)b-xL model class [30] is not 
contained in Eq except for x = 1, and a different anomaly-free completion of the model 
{e.g., involving charged lepton singlets [30]) is needed. A similar remark applies to the 
models in Ref. [12] which predict Z' charges of the SM fermions as in E^ but distinct 
charges of exotics. E.g.., the model [12] couples like the Z^ to the SM fermions as 
well as to the exotic charged leptons, but the D and D charges are multiplied by a 
factor of 3/2 and further SM singlets must be added (see, e.g., Table III in Ref. [3] for 
details and generalizations). 

Finally, we will consider two models corresponding to maximal constructive {Z^_^^^) 
and destructive interference of the Z' amplitude for u quarks (dominating the 

Drell-Yan production process at large momentum transfer) with those of the 7 and the 
ordinary Z boson (see Section 4). 
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Figure 1: The solid red curve corresponds to the CDF 95% C.L. upper limit on the DY 
production cross-section of a Z' boson as a function of its mass. The intersection points with 
the theoretical cross-sections o"^?"^'^^ for various models using Eq. (A. 9) and CTEQ6L PDFs 
give the corresponding Mz' lower bounds. See Table 4 for numerical values obtained by this 
procedure. 



3. Direct searches at hadron colliders 

Due to the large QCD background at the Tevatron, the decay into a lepton pair is the 
preferred discovery channel for a Z' since leptons are relatively easy to identify and their 
energies and momenta can be measured more precisely than those of hadrons, although 
h quarks [33, 34], t quarks [34, 35] and di-jets [36, 37, 38] can also be detected. Among 
leptons, the background for a r pair is harder to manage [39] compared to the 
and e'^e~ channels, of which the former is preferable still [40]. 

The theoretical production cross-section of a Z' boson at a hadron collider depends 
on certain crucial factors, such as the treatment of the parton distribution functions 
(PDFs) of the ingoing quarks and of the radiative corrections to the leading order (LO) 
process. The PDF sets for quarks and gluons are evaluated at various perturbative 
orders for a wide range of factorization scales and momentum fractions by a number of 
independent groups. These sets generally agree with high precision, so that the choice 
depends on whether a particular group provides PDFs at the required perturbative 
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order, the inclusion or neglect of small corrections, the data sets available as of the latest 
update, etc. For the publication [11] we base our analysis on, the CDF Collaboration 
has employed CTEQ6 [41] PDFs. The PDF sets have since been updated a number of 
times by the CTEQ group. We redo this analysis using the latest sets available and 
verify the limits using the MSTW set [42]. 

As shown in Eq. (A.l) of the appendix, for every parton the next-to-leading order 
(NLO) differential Drell-Yan (DY) cross-section consists of three main parts: the PDFs 
for the incoming hadrons, the parton-level hard cross-section and the QCD higher or- 
der terms. The determination of the PDFs requires experimental input. To evaluate 
them, the parameters of some functional form are fit to the data sets from a number of 
experiments (see, e.g., Ref. [43]). The central fit, 5*0, corresponds to the minimum of 
the function. To allow error estimates the CTEQ and MSTW Collaborations also 
provide PDF sets, Sf, which are defined as the eigenvectors of the Hessian matrix [44]. 
Thus, the Sf" are uncorrelated by construction, providing an efficient method of calcu- 
lating the induced variations of the PDF predictions for a chosen practical tolerance 
value, T, defining the region of 'acceptable fits' with Ax^ < T"^- The eigenvectors 
of the Hessian matrix are normalized in such a way that the confidence levels corre- 
spond to hyper-spheres. The uncertainty can then be computed from the simple master 
formula [41], 

1 



AX 



1/2 

(3.1) 



J2[X{St)-X{S, 

where A''^ is the number of eigenvectors, X is the observable (in our case the Z' cross- 
section az') and X{Sf) are the predictions for X based on the PDF sets S^. 

CDF used CTEQ6 PDFs for their calculation of the DY cross-section for which 
CTEQ employed the Particle Data Group average for the strong coupling, as{Mz) = 
0.118 [45]. The CTEQ6M package contains the LO sets in addition to the central NLO 
PDF sets as well as the eigenvector sets for the latter. The latest version 6.6M [46] 
has 22 pairs of eigenvector sets for error calculations. In the PDF sets produced by the 
MSTW group, as has been treated as a fit parameter. This results in LO, NLO, and 
NNLO as{Mz) values of approximately 0.139, 0.120 and 0.117, respectively [42, 47]. 
The MSTW sets contain LO, NLO and NNLO PDFs, in each case along with 20 pairs 
of eigenvector sets. 

The QCD corrections to the LO hadronic process of Z' production may considerably 
alter the magnitude of the cross-section. Conventionally, these corrections are taken 
into account with a '/^-factor', labelled here as i^'"/", and is defined as, 

K"'-'-^ . (3.2) 



dM \ dM ^ 

where cr^'"^'^ is the differential cross-section to order ^ m in a^- This factor expresses 



^ While it may be obvious from the definition (3.2), we recaU that the cross-sections need to be 
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Table 4: 95% C.L. limits on the masses of some benchmark Z' models. Given in the first 
column are the limits which we obtain following the CDF approach as illustrated in Fig. 1. 
We used here CTEQ6L PDFs in order to be able to directly compare our results with those 
published by the CDF Collaboration [11] which are shown in the second column. The third 
column contains the limits obtained from the electroweak precision data [5]. Finally, 95% 
C.L. limits projected for an expected integrated luminosity of 8 fb~^ are listed in the last 
column. 

higher order corrections to Z' production only and not to the complete process. If we 
consider only one quark flavor in Eq. (3.2) then /C™/" is independent of the Z' model 
and Mz'- Thus, the proper way to account for higher order QCD corrections is to 
calculate a different i^"*/" for every flavor, even though it is common practice to choose 
a universal factor for all flavors and models [30] . The main results of the present paper 
are calculated for a^]''^'™^ defined in Eq. (A. 9), with a factor K'^^^ included as we now 
explain. 

As intimated earlier, CDF has used a LO expression (with LO PDFs) and a i^^/o. 
factor (taken from Ref. [30]) for their calculation of the DY cross-section via Z' ex- 
change. We have compared the K^^^ factors from that work with effective K^^^ values 
obtained using the expressions given in the appendix A and found good agreement 
(within a few %). Then we also adopted /('^/^-factors in such a way that our results 
are effectively NNLO. Only for the comparisons in Table 4 with the CDF results we 
use the factor i^^/" instead^. As for the LHC, NLO results suffice [48] and we take 

We take the 95% cross-section upper limit from the data curve in Fig. 3 of Ref. [11] 
and find the intersection with our cross-section, a^™^, as defined in Eq. (A. 9) for 

evaluated with the corresponding order PDFs. 

^For invariant masses beyond 1 TeV at the Tevatron, we use constant factors /v^/'^ = 1.302 and 
if^/i = 1.137. 
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Figure 2: a-/3 Sanson-Flamsteed projection of Z' models. The contours show the 
95% C.L. Mz' hmits in TeV which are from the intersection of o-^?"™^ with the CDF 95% C.L. 
upper hmit on the cross-section. The dotted, dashed, dot-dashed and indigo colored contin- 
uous lines correspond to one-parameter models. For this plot we use CTEQ6L PDFs. 

Z' boson exchange. Final state radiation effects [49] can be ignored in this case since 
we integrate over almost all invariant mass range in such a way that these effects are 
mostly canceled according to the Kinoshita theorem [50]. The evaluation of the multi- 
dimensional integrals involved is done using the CUHRE and SUAVE programs under 
the CUBA package [51]. We input the numerical values of couplings and charges from 
the FORTRAN package GAPP [52]. Only Z' decays into SM fermions are assumed 
and fermion masses neglected. A numerical routine extracted from the PEGASUS 
package [53] is used for the running of a.^ from its value at Mz to the factorization scale 
M. Fig. 1 shows the 'model-lines' for the Z' bosons including some models not included 
in the original GDF analysis. The slopes of the model-lines and their intersection points 
with the experimental data line (giving the Mz' limits) match with those in the original 
CDF plot within a few per mille for models included in both analyses. The 95% C.L. 
mass limits for various models are listed in Table 4. For comparison, the CDF limits 
from [11] and the EWPD limits from [5] are also quoted. The last column in the table 
gives the 95% C.L. mass limits anticipated at the end of the current Tevatron run, 
obtained using the Bayesian statistical method explained in the next section. 
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Z' 


CTEQ6M [41] 


CTEQ6.6M [46] 


MSTW2008 [42] 


Zx 


885 


886 


875 


Zi 


779 


783 


764 




1070 


1068 


1063 



Table 5: 95% C.L. limits on Mz' obtained using three different NLO PDFs. 

Fig. 2 shows the area preserving sinusoidal (Sanson-Flamsteed) projection of a 
hemisphere parameterizing the Z' bosons in terms of the Eq angles a and /3. We note 
that the Mz' lower bound for Z^ji quoted in Table 4 corresponds to the normalization 
given in Ref. [3], and hence differs from the corresponding value in Fig. 2, where all the 
models are normalized as in Eq. (2.4). 

In Table 5 we display the dependence of these limits on the PDFs. To also in- 
vestigate the uncertainties due to them we use the central CTEQ6M PDFs and the 
corresponding eigenvector PDF sets as displayed in Table 6 for two selected values of 
Mz'- One sees that the relative uncertainty in az' is very large for the Zj which is due 
to Q'^ = Q'^ = in this model and points to the fact that for a given M the uncer- 
tainties in the d quark PDFs are larger than those of the u quarks. Also, the d quark 
contribution is suppressed by more than an order of magnitude with respect to that of 
the u quark which is reflected by the ratio of the cross-sections for the Z/ and the Z^. 
We recall that we use a common normalization for all models so that the cross-sections 
can be directly compared. Finally we show in the table how the uncertainty in az' 
affects the Mz' limits in these models. 

4. The Bayesian statistical method 

The CDF Collaboration collected an integrated 2.3 fb~^ of data [11] in the fi'^fi~ 
channel, binned in inverse invariant di-muon mass, m'^. The CDF analysis then looks 
for an enhancement in di-muon production above the SM background for particular 
Eq models, and so their lower limits on Mz' correspond to upper limits on the cross- 
section^. But Z' bosons interfere with the SM neutral gauge bosons, and destructive 
interference would result in a reduction of the SM cross-section. For this reason in 
addition to the general motivation given in Section 1 — clear-cut combination of Z' 
constraints from quite distinct sources — we adopt a statistical framework wherein it 
is straightforward to address interference effects and to vary the coupling strength (see 
also Ref. [54]) up to the strong coupling regime (and broad resonances). 

The basic idea is to apply the Bayesian analysis of the SM Higgs mass, Mh, of 
Ref. [55] to Z' physics. In this case the collider constraints from LEP 2 [56] and the 

^It utilizes signal templates that have been generated with a fixed and relatively narrow Tz' = 
2.8% X Mz' (motivated by Tz). 
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Mz' = 0.8 TeV 


Mz' = 1.1 TeV 




Z' 


fTz' [fb] 


Aa 


cTz' [fb] 


Aa 


Mz' [GeV] 


Zi 

Z4 


8.6 ±0.7 
3.0 ±0.7 

54 ±4 


8% 
24% 
8% 


0.41 ±0.04 
0.14 ±0.03 
2.6 ±0.3 


9% 
20% 
11% 


885 ±9 
779^1 

1070+}^ 



Table 6: Di-muon cross-section for Z' exchange obtained using NLO CTEQ6M PDFs sets 
for two representative values of Mz', along with the uncertainty due to the eigenvector PDF 
sets. The percentages of the uncertainties are given in the third and fifth columns. In the 
last column we give the 95% C.L. limits on Mz' and the positive and negative uncertainties 
in these. We obtain the uncertainties from the intersection of a'^'^'^^^ ± Ait^V*^"^^^ for Z' 
exchange and using Eq. (3.1) with the CDF 95% C.L. upper limit on the cross-section. 



Tevatron [57] were included using the published log-likelihood ratios, 

, ^ , io(data|s±6) , 

LLR..-21„£,.-21„«-|j-L_i. (4.1) 

These are given in terms of the probabilities (likelihoods) to obtain the data, condi- 
tional on the signal plus background hypothesis, p(data|s ± b), and background only 
hypothesis, p(data|6), and may be compounded of many experiments, channels, ener- 
gies, etc. The £j depend on the parameter(s), /i, of interest, (/i = Mh in Ref. [55]), 
through the signal hypothesis. Information on fi is obtained by Bayes's theorem, 

,(,a..a)^l^£^, (4.) 

where p(/i) is the prior probability density function (pdf) entering Eq. (1.1), and is a 
summary of our knowledge, if any, prior to the experiment or analysis. In the absence 
of prior information, or if the prior information is explicitly taken into account by extra 
factors p{fi\ni), then p{fi) is called non-informative, and is most conservatively taken 
as p(/i) = 1 or p{fi) = fi~^, whenever fi may be an arbitrary real number or positive 
real number, respectively. Notice, that j9(data) drops out from likelihood ratios. This 
is crucial: if various data points show poor compatibility, or if an excess (or deficit) 
is observed, this will have an impact only if some value of fi describes the data better 
than some other. 

We use here a similar philosophy as in Ref. [55], but here we first have to construct 
the corresponding LLR ourselves. Our input data are the number of events, Ui, in 
bin i (see Fig. 1 in Ref. [11]). The parameter set /i may include all of the parameters 
introduced in Sections 1 and 2, namely Mz', g', Ozz', a, and (3, but in this paper we 
allow only Mz' and g' for some specific models {i.e., fixed values of a and /3) and set 
0ZZ' = 0. Complementary data sets, such as other channels, LHC and D0 results, 
EWPD data [5, 6], and LEP 2 constraints [7], will be necessary to disentangle these 
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Figure 3: The x dependence of the (proton) PDFs at = = (1000 GeV)^. The upper 
(lower) plot shows the LO (NLO) PDFs of u, d, u, and d quarks provided by the CTEQ and 
MSTW groups. 

parameters in an integrated analysis. The SM point corresponds to M^/ = or g' = 0. 
The events, Ui, in each bin follow Poisson statistics, 

Pin,\u,) = (4.3) 
nil 

where Ui is the predicted number of events in bin i given specific values for Mz' and g' . 
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It is important to note here that the above hkehhood is the same as the one 
employed by the CDF collaboration in their analysis [58]. However, they determine v'^ 
by summing the bin-counts from the SM and the bin-counts from the signal template 
as mentioned earlier, effectively summing the cross-sections of the SM process and the 
Z'-mediated process without any interference. Therefore, the CDF approach essentially 
differs from ours not in the choice of the likelihood, but in that the interference effects 
were neglected in their analysis in order to keep it simple and model-independent [59] , 
deliberately making it blind to wider resonances through the use of templates based on 
narrow signal width. We, conversely, treat coupling strength as a free parameter and 
avoid signal templates, which makes the inclusion of interference effects rather natural 
in our framework. 

In practice, we compute a grid^ of values for the ui discretizing Mz' and g' and 
interpolate in every one of the first 35 bins corresponding to the invariant mass range 
searched by CDF. We have thus effectively reduced the analysis to a least-x^ fit where 
any observed event count adds a piece. 

Ax? = LLRi = 2(^[-Ui + m In , (4.4) 

to the overall function (i/j and v[ refer here to the SM and SM plus new physics 
expectations, respectively). 

The detector resolution, A = 0.17 TeV~^, is approximately constant (in the vari- 
able ?Ti~^), and must be taken into account since it is of the order of the bin size of 
(3.5 TeV)~^. We define M^^ as the theoretical invariant mass of the muon pairs, as 
opposed to the nominally measured m^^, and introduce the convolution. 



oo 

C dm-lA{m,^) j dM;^p{m-l\M;^)K^/^K^^^^+Unu. 



(4.5) 



where C = 2.3 fb"'^ is the integrated luminosity, e is the detector efficiency, which we 
take as a constant 0.982 for all bins, a^^^ is given by Eq. (A.l), and Unm- refers to the 
non-DY background which is extracted directly from Fig. 1 of Ref. [11]. A{m^^) is the 
total acceptance of the CDF detector, which increases from about 0.13 at the Z-pole 
to about 0.4 at 1 TeV and then falls off rapidly [11]. For our current analysis, the 
acceptance values have been gleaned from [60]. 

Along with the QCD corrections, the 0{a) QED radiative corrections [49] also 
have a sizable effect on the DY cross-section, and strongly affect the shape of the 



^Alternatively calculating directly without the grid gives mass limits which differ by at most 
3 GeV for our benchmark models. However, this dramatically increases CPU time if a multi-variate 
minimization is performed i.e., without fixing g' and other model parameters. It is therefore expedient 
to avoid the mostly redundant PDF integrations. 
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Figure 4: The 95% C.L. exclusion contours in the g' vs. M^/ plane for the Z^_-^^ and 
bosons. These are compared to the contours that are obtained when the interference effects 
are ignored. 

di-lepton invariant mass distribution. While initial state radiation is negligible for di- 
muon masses between 50 and 100 GeV, final state QED corrections are in fact larger 
than the 0{as) QCD corrections and so have to be taken into account. The QED 
corrections are shown in Fig. 6 of Ref. [49] with a rapid variation visible in the range 
40 GeV < M < 110 GeV. Just below the Z peak, these corrections enhance the 
cross-section by up to a factor of 1.9, so the cross-sections in the neighboring bins 
to the Z peak differ considerably from the values expected without these corrections. 
When the full di-muon mass range is integrated over, the large negative and positive 
corrections tend to cancel and do not have a big impact on the total cross-section [50] . 
For M > 130 GeV, they uniformly reduce the differential cross-section by 7%, and for 
our calculations^ we have extrapolated the data points from the mentioned figure and 
used these as multiplicative factors, which we refer to as K^{M^^). In principle, such 
effects should also appear near the Z'-po\e, but considering that the bin around the 
expected Z' mass is fairly wide, any large effect around the peak will be washed out at 

^The figure [49] has been generated for i/s = 1.8 TeV, while our process is being computed for the 
^/s — 1.96 TeV of the current Tevatron run, but we expect this to have negligible effect on our final 
results. 
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Figure 5: Contours in the g' vs. M^, plane for the model. The solid lines correspond — 
from top to bottom — to the 95%, 85%, and 68% C.L. contours using the Bayesian analysis, 
and are to be compared with the dotted line using the CDF (frequentist) approach. The 
dot-dashed line shows what one would expect from the Bayesian method based on pseudo- 
experiments (PE), while the EWPD yield the dashed line. The horizontal line indicates the 
value g' = 0.461 which is motivated by gauge coupling unification and often used as reference. 
The cross indicates the best fit value. 



least for weak and intermediate coupling strength. 

Returning now to Eq. (4.5), the quantity p('ti~^|M~^^) is our smearing function, 

Pi'nfi-lW;^) = m^^ (4.6) 

with = M^^A^ and = m^^j,M^^j, and is constructed as a Beta distribution 
with mean M^^ and variance A^. Note that it approaches a Gaussian form for a ^ 1 
(which is the case except for the first few bins), but is more adequate than a Gaussian 
since M^^ takes non-negative values only. We note in passing that we neglect here sys- 
tematic and theoretical uncertainties, justifying this with the very small event numbers 
in the most relevant bins so that statistics dominates. 

Eq. (4.4) also makes it explicit how our approach allows the new physics to enter 
with either sign, as is always the case for interfering Z' bosons. At the level of the 
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differential cross-section, tlie interference terms cliange sign wlien M^^ crosses tlie Z 
or Z' poles. Thus, there are fairly large cancellations at work when the whole range 
of m^^ is integrated over, and when the objective is the usual hunting for a narrow 
bump where neglect of interference effects is justified [61]. Since here we put more 
emphasis on the event distribution over larger numbers of bins, the interference issue 
becomes more interesting. To have a closer look as to how significant the interference 
effects are numerically, we now discuss two cases where they are enhanced, and 
They are defined to have, respectively, maximum constructive and destructive 
interference with the SM amplitudes for up quarks in the limit M^^ — )■ oo, i.e., we 
extremize the expression (in a slightly more compact notation), 

e'QuQAQLu + Q'r.Ml, + Qr,) + 91{^LuQ'lu + ^RuQru){^l^Q'l, + ^R^Qr^)- (4.7) 

We chose = as our reference value because then any dependence on Mz' and 
V z' drops out. Moreover, at large M^^ the down quark contribution to the PDFs is 
strongly suppressed, providing a further simplification. For the case of the it now 

turns out that neglecting the second (ZZ' interference) term shifts the corresponding 
values for our model parameters a and /3 only at the 10~^ level, so that we can neglect 
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this term, as well, and we find. 



1 

a = - arctan 
2 



2^6 



71 



76 



■ o 



/3 = 



(4.8) 



9 



2 



which is relatively close to the Zb-l case. The facts that the Zb-l boson couples 
only vector-like and that the vector-coupling for the muons is suppressed by a factor 
1 — 4sin^ 9\y{Mz') ~ 0.04, may give a rationale for why in this case the ZZ' interference 
term is small. Similarly, the Z~_-^^ is numerically close to the Z^ boson which has only 
axial-vector couplings to the SM fermions. In this case we can neglect the first {'fZ' 
interference) term in Eq. (4.7) and simply define Z~_-^^ = Z^. As for the integrated 
cross-sections, the constructive interference for Z^_-^^ is about an order of magnitude 
larger than the destructive interference in Z~_^^^, and in the latter case we find that 
the sign of the interference effect in the total cross-section is reversed compared to the 
amplitude level in the M^^j — )■ oo limit. 

We illustrate the interference effects for the Z^_^^^ and Z^ bosons in Fig. 4. As 
can be seen, they become significant for g' values of order unity. In fact, for large 
Mz' ^ 0{1 TeV) most of the expected signal events come from the 7Z' interference, 
since the pure Z' exchange contribution is more strongly mass suppressed. Another 
way to quantify the interference effects is to look at the behavior of the best fit location. 
E.g., for the Z^ model (not included in the plot) we found the global best fit at Mz' = 
2.05 TeV and g' = 1.14, with the value of Ax^ = -1.43 relative to the SM. On the 
other hand, if we turn off the interference effects, the global minimum strongly shifts 
to Mz' = 0.189 TeV and g' = 0.011 with Ax^ = -0.97. 

We also stress that the interferences are important if one wants to discriminate 
between models^. E.g., the Z^_-^^ and Z^ models have their global minimum at low 
mass and weak coupling similar to the values above regardless of the interference. But 
the x^-minimum becomes deeper in the presence of interference effects, even though we 
show in Fig. 4 that the mass limits are unaffected at small coupling. Moreover, without 
interference is virtually degenerate at the minimum for the three mentioned models, 
but this is lifted by the interference effects. 

The contours in the g' vs. M'^} plane are given in Fig. 5 for the Z^. For comparison 
with our approach, we extended the CDF limit corresponding to g' = 0.461 to other 
values of the coupling. Crucially, the CDF line breaks down at g' ^ 1.5 owing, again, 
to the fact that their templates assume a narrow Z' boson. As can be seen, our method 
reveals a strong variation of the Z' mass limit with g', while the frequentist method used 
by CDF shows a weaker and mostly monotonic dependence. The PE line is obtained 
by assuming an event count rii = Ui {i.e., the SM expectation even though Ui is not 
integer valued) in Eq. (4.4), in place of the actually observed number of events for all 
the bins. This yields a smooth and monotonic contour, demonstrating that the strong 

^The importance of interference effects in forward-backward asymmetries was emphasized in 



Ref. [62]. 
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Z' 


two-parameter 


one-parameter 


frequentist 


z^ 

A. 


886 


914 


885 




888 


916 


871 




913 


939 


898 


Zi 


784 


809 


779 


Z(Ji 


1043 


1064 


1070 


7+ 

u—int 


1056 


1084 


1110 



Table 7: 95% C.L. limits on Mz' [in GeV] for six representative Z' bosons. Shown are the 
two-parameter (with g' free) and one-parameter (with g' = 0.461 fixed) hmits. In the last 
column we show the one-parameter limits obtained using the frequentist method as discussed 
in Section 3. 

Mz' dependence of the g' limits is real. We also show the contour g' oc Mz' from the 
EWPD. As exemplified here, the EWPD give stronger (weaker) constraints for larger 
(smaller) values of g' when compared to the CDF data. 

At a crude level, the exclusion curves obtained by us using either the CDF method 
or our own are quite similar, and one notices that the downward fluctuations in the 
cross section data [11] at around 3 and 6 TeV~^ are reflected in both cases. At a 
flner level, we tend to see slightly higher Z' mass limits using our Bayesian approach 
for small g' , compared to those obtained using the CDF approach, while for stronger 
coupling, g' > 0.5, the limits are lower for most models. Our method also maps out 
fluctuations more faithfully. These small differences arise from a combination of effects, 
such as the interference term, the statistical interpretation, and to a smaller extent our 
neglect of systematics. The best flt is also indicated in the flgure. It is amusing that it 
occurs for Mz' values close to but somewhat smaller than the invariant masses where 
the CDF Collaboration [63] sees a signiflcant deflcit (at rngg ~ 200 GeV) followed by a 
signiflcant excess (at ~ 240 GeV) in di-electron DY production. 

As mentioned earlier, it is not possible to extend the CDF limit to values of g' 
much greater than unity. Thus, our approach provides a convenient way to flnd mass 
limits for a strongly coupled Z' boson. The Mz' limits for a particular g' value can be 
read off from Figures 4, 5 and 6. Note, however, that these are two-parameter limits 
since there are two flt parameters in the minimization of the function. The 95% 
C.L. in this case corresponds to Ax^ = 5.99, while a one-parameter limit with g' flxed 
corresponds to Ax^ = 3.84. To illustrate this we compare in Tab. 7 the values of the 
two-parameter and one-parameter limits from the Bayesian method corresponding to 
g' = 0.461, and also show the limits of the (one-parameter) frequentist method for 
various models. As can be seen, the variation of the limits due to different statistical 
interpretations is larger than any of the theoretical uncertainties reviewed in Section 3. 

Finally, we exploit the predictive nature of the Bayesian formalism to give a general 
idea of how various collider options compare. For this we project 95% C.L. lower limits 
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2 TeV 


7 TeV 


14 TeV 


28 


TeV 


/£[fb-i] 


pp 


pp 


pp 


pp 


pp 


pp 


pp 


pp 


3 


0.65 


0.88 


1.60 


2.11 


2.46 


3.26 


3.76 


4.52 


10 


0.74 


0.98 


1.86 


2.50 


3.07 


4.05 


4.58 


5.87 


30 


0.82 


1.08 


2.14 


2.85 


3.56 


4.52 


5.36 


7.15 


100 


0.90 


1.20 


2.43 


3.22 


4.10 


5.29 


6.23 


8.59 


300 


0.97 


1.25 


2.69 


3.55 


4.60 


5.97 


7.85 


9.89 


1000 


1.05 


1.32 


2.98 


3.89 


5.14 


6.70 


8.98 


11.2 


3000 


1.12 


1.38 


3.24 


4.18 


5.63 


7.33 


10.0 


12.5 



Table 8: Projected 95% C.L. exclusion limits [in TeV] on Mz' for the model using our 
Bayesian method. These limits are obtained by assuming the number of observed events rij 
equal to the SM expectation Ui in Eq. (4.4). We consider typical CM energies (shown in the 
top line) and a range of integrated luminosities / £ for pp and pp colliders. See the text for 
details. 

on the mass of the boson for various reference CM energies and luminosities (actual 
and hypothetical) for pp as well as pp collisions. We employed the CDF acceptance 
(asymptotically for large dimuon masses, we used the constant value of 0.316) and 
ignored FSR as well as finite resolution effects in the dimuon invariant mass. For the 
resulting limits listed in Table 8 we fixed the number of bins to ten, and varied the bin 
size until the mass limit reaches a maximum. This optimal bin size turned out to be a 
constant times the inverse of the corresponding limit. The cases of 3 fb~^ and 10 fb~^ 
correspond roughly to the currently analyzed and final Tevatron data sets, while 30 fb~^ 
refers to what is often called the low-luminosity LHC. Likewise, 300 fb~^ and 3000 fb~^ 
correspond to the high-luminosity and luminosity-upgraded LHC, respectively. The 
cases of 100 fb~^ and 1000 fb~^ are included so that they can be compared with the 
existing projections in the literature [64, 65, 66]. The limits in Table 8 are much higher 
for pp colliders where the q mostly emerges as a valence quark from the anti-proton 
with a PDF given by = while in the case of pp collisions the q is always a sea 
quark and thus has lower PDF values (as illustrated in Fig. 3). Table 8 shows that for 
Z' (DY) physics pp colliders have a relative advantage as significantly less integrated 
luminosity is needed to match the pp case, although larger CM energies tend to mitigate 
this effect. 

5. Conclusions and outlook 

In this article we have completed a step towards an integrated analysis of the physics 
parameters associated with extra Z' bosons. These parameters, Mz', g', Ozz', (3, 
etc., need to be disentangled and this can be achieved by exploiting the complemen- 
tary nature of electroweak precision data, lepton and hadron colliders, as well as the 
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discriminatory power of simultaneously analyzing a variety of processes within each of 
these data classes. The (a, (3) parameterization [67] discussed in Section 2 allows us to 
analyze fairly different kinds of models, such as Eq derived models with and without 
kinetic mixing, Z' bosons motivated by chiral models of weak scale supersymmetry [12] , 
models based on left-right symmetry, and other one-parameter models, all on the same 
footing. 

To move beyond a collection of lower mass limits on the Z' bosons and towards 
an integrated analysis it is crucial to have a common framework. We proposed such a 
framework in Section 4 after having laid out some technical groundwork in Section 3, 
including a detailed formulation of the DY production of fi~^fi~ via neutral gauge bosons 
at hadron colliders such as the Tevatron and the LHC, and a discussion of the asso- 
ciated uncertainties. We have shown that the limits obtained with our approach are 
numerically close to those in the more traditional frequentist approach. This should not 
come as a surprise. Indeed, while the need to introduce a prior probability density is 
often held against Bayesian data analysis (a point which has been over-emphasized by 
many authors), the posterior distributions are often (and certainly here) strongly data 
driven. Our resort to Bayes originates more out of the general philosophical mindset it 
represents, and the ideal application it provides to problems of parameter estimation. 
A welcome by-product of this approach is the ease with which the U{iy coupling can be 
varied, the interference effects can be included, and the exclusion limits for anticipated 
future runs of various collider experiments can be predicted. 

We are looking forward to carry out similar analyses with the first results by the 
CMS [9] and ATLAS [10] Collaborations. There will also be further results (see, e.g., 
Refs. [68, 69]) and different channels by the CDF and D0 Collaborations. Finally, a 
more systematic classification and discussion of the model class of Section 2 is also 
underway [19]. 
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A. Cross-sections and PDFs 

The NLO differential cross-section for the DY process with a neutral gauge boson G 



- 21 - 



as the mediator, pp — )■ GX — )• /+/ X, is given as [70], 



da 



NLO 



dM 



M / dzdxi e 1 

NcS J xiz \ 



Xizrl 



(A.l) 



+ /f(x2,M^)[/^^(xi,M^) + /^^(xi,M^)][ X 



27r 



where A^c = 3 is the color factor, M is the invariant mass of the observed lepton pair 
and y/s is the energy of the pp colhsion in the CM frame, = y/s/M, and ^ = xizr^. 
f^l^ are the PDFs of the quarks and gluons coming from hadron A. is the strong 
couphng constant, and 
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with the Cp = 4/3 and Tr = 1/2, and the '+' distribution defined as 



dzg{z] 



log(l - z) 

1 - z 



dz{g{z)-g{l)}{^^ 



logfl - z) 



(A.3) 



The expression for the hard scattering cross-section of the process qq — ?■ i^i 

da 



aqq^e+i-{M 



I dcos9* 



d cos 9* 



(A.4) 



= {\All\' + \Annn (1 + cos 9^ + (I^.^P + I^^lP) (1 - cos 9^}, 

is given in terms of the polar angle, 9*, in the CM frame and the individual amplitudes. 



Aj = -Q{q)e^ 



gle,iq)e,{i)M^ 



g^Q'MQ'MM' 



- M} + iMzTz - Mj, +iMz'Tz' 



(A.5) 



where i,j run over L,R. Q{q) is the electric charge of the quark and e = gsin9\Y- 
Mz,z' and Tz,z' are the masses and total decay widths of the Z and Z' bosons. 



6l(/) = Tsif) - Qif) sin^ 9w, enif) = -Q{f) sin' 9 



w, 



(A.6) 



are the effective couplings of the ordinary Z to fermion / entering with coupling 
strength, gz = g/ cos^pi^ = 0.7433. The Q'^if) are given in Table 3 and Q'^i^f) = 
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—Q'l^if). As for the Z' coupling strength, one often employs the (one-loop) unification 
value [16], g' = ^/5/3 sin O^gz = 0.4615, but we will do so only when comparing results 
to other analyses, since one of our goals is to study the g' dependence. 
Integrating Eq. (A.l) gives the LO differential cross-section^ as 



da 



LO 



dM 



■-—M / dxi 



Xi 



a{M'){f^{x„M')f^^ix„M') + f^{x„M')f^ix„M')], (A.7) 



where X2 



now. Note here that for purely SM contribution to the process, the 



third term, while for the expected cross-section via a Z' boson only, the first and second 
terms in Eq. (A. 5) can simply be ignored. The decay width, Tz', given in eq. (A. 5), is 
the sum of the partial decay widths of the Z' boson into all the fermions it couples to. 
The partial decay width into a Dirac fermion pair is written as [71] 



g"Mz 
2An 



1- 



4M| 
Ml 



Ml 
Ml 



(Ql + Q'l) 



6M] , / 



M2 



where Mf is the mass of the final-state fermion. We add the factor M"^ /M"^, to get an 
's'-dependent Z'-width [72]. For the range of Mz' of interest here, Mf -C M^/ for SM 
fermions, and the above expression becomes independent of the fermion masses. 

With the above definitions we always obtain an effective NNLO result by multi- 
plying by the corresponding K"^^^ factor in Eq. (3.2). For NLO and LO PDFs the 
corresponding effective NNLO cross-sections are, respectively. 



„NLO-PDF 

az> 



dMK^/ 



I da 



NLO 



dM 



„LO-PDF 

azi 



(A.9) 



where M„-jin = 0.1 TeV is the lower invariant mass of the analyzed events in [11]. 
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